gr-qc/9708036 
BUTP-97/23 
August 17, 1997 



Non-Abelian black holes: 
The inside story ^ 

Peter Breitenlohner George Lavrelashvili ^ Q 
and Dieter Maison § 

Max-Planck-Institut fiir Physik § 
— Werner Heisenberg Institut — 
Fdhringer Ring 6 
80805 Munich (Fed. Rep. Germany) 

Institute for Theoretical Physics ^ 
University of Bern 
Sidlerstrasse 5 
CH-3012 Bern, Switzerland 

Abstract 

Recent progress in understanding of the internal structure of non-Abelian 
black holes is discussed. 



^Talk given at the international Workshop on The Internal Structure of Black Holes 
and Spacetime Singularities, Haifa, Israel, June 29 - July 3, 1997 

^On leave of absence from Tbilisi Mathematical Institute, 380093 Tbilisi, Georgia 



1 Introduction 



Let us start recalling what is known about the non-Abelian black holes in 
general. The non-Abelian story begun in 1988 when Bartnik and McKinnon 
(BK) unexpectedly found flj numerically a discrete sequence of globally reg- 
ular, particle like solutions of the Einstein- Yang-Mills (EYM) theory. Soon 
the same model was solved with the different boundary conditions corre- 
sponding to black holes 0. Numerical findings were confirmed by math- 
ematically rigorous existence proofs 0, Q of both regular and black hole 
solutions. It turned out that all the above solutions of the EYM theory are 
classically unstable against small perturbations. In addition to the genuine 
gravitational instabilities there also instabilities of topological origin p, ^ 
related to the sphaleron nature of the solitons [||, H. 

A few related systems were investigated. It was shown 



10, 11 that the 



gravitational field can be replaced by a dilaton, so that the YM-dilaton the- 
ory in flat space has a tower of solutions similar to the BK sequence. The 
combined EYMD theory |l^, [l^, [l^, [l^, [T^ was shown to possess both reg- 
ular and black hole solutions for any value of the dilaton coupling constant. 
For the EYMH theory with a Higgs doublet it was shown that the the- 
ory in addition to the gravitating sphaleron solution contains its BK type 
excitations. 

Furthermore the EYMH theory with a triplet Higgs [[Tsj |19|1 was studied. 
This theory is interesting since in the flat limit it contains t'Hooft-Polyakov 
monopoles, which are known to be stable. It was shown that the basic 
monopoles continue to exist, when gravity is switched on, at least as long as 
the gravitational self-interaction is not too strong. In addition the monopole 
admits unstable BK type excitations. 

Adding a cosmological constant to the EYM theory one obtains non- 
asymptotically fiat analogues of the BK solutions |]2D|, pT| . 

Whereas in the above study the globally regular solutions were completely 
analyzed, the investigation of the black hole solutions was not complete since 
their internal structure was unknown. Recently this problem was investigated 
independently by us ||2^ and by Donets, Gal'tsov and Zotov ||22|. Our main 
results on the EYM case essentially agree with theirs, although we differ in 
some details. 

We found it adequate to describe the generic behavior as a kind of mass 
inflation closely related to the "usual" mass inflation (see e.g. p4|, |25|, ^ 
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28| and many references in the present proceedings). 

In addition to the generic behaviour there are three different types of 
special solutions, which are obtained by fine tuning of the initial data at the 
horizon. There are solutions with Reissner-Nordstr0m (RN), Schwarschild 
and pseudo-RN type behaviour. 

The present contribution is essentially based on our paper p3|. 



2 Field Equations 

The action of the EYMH theory is 




where g denote the gauge coupling constant, G is Newton's constant, F is 
the field strenght of the SU(2) Yang-Mills field and V{^) is the usual quartic 
Higgs potential. The pure EYM action and corresponding equations can be 
trivially obtained from the EYMH ones by putting the Higgs field $ and its 
potential V{^) to zero. 

For the static, spherically symmetric metric we use the parametrization 

ds^ = A^Bdf - ^ - r\R)dn'' , (2) 

with dVL^ = dO"^ + siv? 6d({P' and three independent functions A, 5, r of a 
radial coordinate which has, in contrast to r, no geometrical significance. 
As long as dr/dR ^ the simplest choice for r is = r, i.e. Schwarzschild 
(S) coordinates. In this case it is common to express B through the "mass 
function" m defined by i? = 1 — 2m/ r. 

For the SU(2) Yang- Mills field VF^ we use the standard minimal spheri- 
cally symmetric (purely 'magnetic') ansatz 

ly^T^dx'^ = W{R) {Tide + T2 sin Odip) + T3 cos Od^ , (3) 

and for the Higgs (triplet) field we assume the form 

•l-X = H{R)n''Ta , (4) 

where Tq denote the generators of SU{2) in the adjoint representation. One 
might also consider other representations for the Higgs , e.g. doublet, but we 
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believe that the behaviour near r = is the same. Plugging these ansatze 
into the EYMH action results in 



with 



S = - / dRA[l{l + B([rr- + ^-^(.r')')) - BrV, - Vi] , (5) 



y^.i^^hnr, (6) 



r2 ■ 2 

and 

2^2 



V, = ^^-^ + ^{H' - a'f + W'H' . (7) 

Through a suitable rescaling we have achieved that the action depends only 
on the dimensionless parameters a and j3 representing the mass ratios a = 
MwVG/g = Mw/gMpi and p = Mh/Mw {Mh and Mw denoting the Higgs 
resp. gauge boson mass). 

Using S coordinates the field equations obtained from (|]) are 

{BWy = W{—^^ — + H^) - 2rBW'Vi , (8a) 

{r'^BHy = {2W^ + _ a^))H - 2r^BH'V, , (8b) 

(rB)' = 1 - 2r^BVi - 2V2 , (8c) 
A' = 2rViA . (8d) 

If dr/dR = (equator!) S coordinates become singular and one has to use a 
different choice (gauge) for R. A convenient possibility is given by i? = 
for B > resp. B = — for B < 0. We denote this radial coordinate by 
r in order to distinguish it from the S coordinate r. With this choice the 
metric takes the form 

ds' = -^dt^ + r\T)idT^ - dn^) . (9) 
r^(r) 

We will refer to this system of coordinates as isotropic coordinates. The 
equations obtained with isotropic coordinates are given in p3[| . 
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3 Singular points 



Obviously the field Eqs. (§) are singular at r = 0, r = oo and for points where 
B vanishes. Note that a generic solution should have five free parameters 
for the EYMH case resp. three for the EYM case. Here we have ignored the 
variable A, which decouples (compare Eq. (|]d)). 

At the singular points the generic solution is singular and only a non- 
generic subset of solutions stays regular. 

In the vicinity of r = cx3 one finds a 3-parameter family of asymptotically 
flat solutions for the EYMH theory (2-parameter family for the EYM case). 

At r = r/j for any given > one finds a 3-parameter family charac- 
terized by the value of a gauge and Higgs fields at the horizon, Wh and 
respectively. 

In case of the EYM theory it was shown p, 0, Q that for any given > 
there is a discrete set of solutions interpolating between the horizon and 
infinity. They can be characterized by an integer n, the number of nodes of 
the gauge amplitude W. 

The EYMH case is more complicated [0, |Ty] . Let us here only recall that 



the previous research was concentrated on the investigation in the interval 
[rh, oo]. 

Now we turn to the classification of the singular behavior at r = 0, which 
is of particular relevance for the internal structure of black hole solutions. 
We have to distinguish two cases, B > and B < 0. 

1. -B > 0: For black holes this case is only possible, if there is a second, 
inner horizon. One finds a 5-parameter, i.e. generic, family of solutions. 



^(^) = ^0 + TT7T-^^1F2T^' + ^^r' + O(r^) , (10a) 



H{r)=Ho + Hir + 0{r') , (10b) 

B(,).M^-?^^ + o(i). (100 



r 



where Wq ^ 1. 

According to the asymptotics of i?(r) we may call the singular behavior 



to be of RN-type. The special case W"^ = l,Mo < leads to S-type 



behavior with a naked singularity. On the other hand Wq = 1, if( 
0, Mo = gives regular solutions. 
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2. B < 0: This case is more involved, with two disjoint famihes of 
singular solutions. 



2.1 There is a 3-parameter family of solutions with a S-type singularity, 
characterized by the asymptotics 

W{r) = l + W2r^ + 0{r^) , (11a) 
H{r) = Ho + 0{r) , (lib) 

5(r) = -^ + 0(l) . (11c) 
r 

where Mq > 0. 

Obviously the condition B < {Mq > 0) prevents the existence of 
regular solutions in this case. 

2.2 There is an additional 2-parameter family of solutions with a pseudo- 
RN singularity (pseudo because -B < 0). 

W{r) = Wo±r + 0{r^) , (12a) 
H{r) = Ho + 0{r^) , (12b) 

(Wn - 1)^ AWnd - Wn) 

B{r) = - ^ ^ ' ±^11^ ^ + 0(1). (12c) 

with Wq 7^ 1. The eigenvalues of the linearized equations are Ai,2 = 
-1/2(3 ±iv^5) and A3 = -1. This is a repulsive focal point which will 
turn out to be important for the cyclic behavior in the EYM case. 

Note that the corresponding Taylor series for the singular points in case 



of the EYM theory were listed (with minor mistakes) in |2^, whereas in our 



work ||23| the more general EYMH theory was studied and local existence 



proof was given. In other words we have shown that the expressions above 
are in fact the beginning of a convergent Taylor series for H and r^B. 

In the case i? < we obtained no singular class that has enough param- 
eters (three for EYM and five for EYMH) to describe the generic behavior. 
Since, also the appearance of a second, inner horizon is a non-generic phe- 
nomenon, one may wonder, what the generic behavior inside the horizon near 
r = looks like. This situation is shown schematically on the Fig. |I] for the 
case of the EYM theory. 
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4 Numerical results 



In order to investigate the generic behavior of non-Abehan black holes inside 
the event horizon, we integrate the field Eqs. (|^) | from the horizon assuming 
-B < 0, ignoring the constraints on the initial data at the horizon required 
for asymptotic flatness. 

Using the Killing time t the horizon is a singular point of the equations. 
Consequently one has to desingularize the equations in order to be able to 
start the integration right there. How this can be done, was described in 

MM- 

As one performs the numerical integration one quickly runs into problems 
due to the occurrence of a quasi- singularity, initiated by a sudden steep raise 
of W and subsequent exponential growth of B resp. m (compare Figs. ||, 
and 1^ for some examples). This inflationary behaviour of the mass func- 
tion is similar to the one observed for perturbations of the Abelian black 
hole solution at the Cauchy horizon ^ |2^. While this fast growth 
continues indeflnitely for the EYMH system, it comes to a stop without the 
Higgs fleld. The mass function reaches a plateau and stays constant for a 
while until it starts to decrease again. When B has become small enough, i.e. 
the solution comes close to an inner horizon, the same inflationary process 
repeats itself. Generically this second "explosion" is so violent (we will give 
estimates on the increase of m in chapter |^) that the numerical integration 
procedure breaks down. 

Besides these generic solutions there are certain families of special solu- 
tions obtained through suitable flne-tuning of the initial data at the horizon. 
There are two classes of such special solutions. The flrst class are black holes 
with a second, inner horizon, the second are solutions with one of the sin- 
gular behaviours at r = for B < described in chapter 0. The numerical 
construction of such solutions is complicated by the fact that both boundary 
points are singular points of the equations. The strategies employed to solve 
such problems are well described in the paper on gravitating monopoles [|T^]. 
Actually, in order to control the numerical uncertainties we used two differ- 
ent methods, which may be called "matching" and "shooting and aiming". 
For matching we integrate independently from both boundary points with 



^ At least part of the solutions possess an equator, i.e. a local maximum of r. For 
those the use of S coordinates is excluded and we integrated the equations in the isotropic 
coordinates. 
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regular initial data, tuning these data at both ends until the two branches 
of the solution match. For shooting and aiming we integrate only from one 
end and try to suppress the singular part of the solution at the other end by 
suitably tuning the initial data at the starting point. 

Our results concerning special solutions are shown in Fig. |^ and dis- 
cussed in p3[. As already said, the first class of special solutions consists of 



black holes with a second, inner horizon; we call them non-Abelian RN-type 
(NARN) solutions. We have determined two such 1-parameter families for 
the EYM system, shown in Fig. |^. As may be inferred from Fig. ^, the (dot- 
ted) curve 2 corresponding to one such family intersects all (solid) curves 
describing asymptotically flat solutions except the one for n = 1. Curve 
continues straight through the parabola rt = I — and runs all the way 
to Th = Wh = 0. The branch to the left of the parabola cannot be obtained 
using S coordinates since solutions develop local maximum of r between the 
horizons. 

Our second NARN family (curve 5 of Fig. ^ stays completely to the left 
of the parabola and ends at rt ~ 0.9 close to the curve 3, whose significance 
will be explained below. 

The second class are solutions without a second horizon (i.e. B stays neg- 
ative) approaching the center r = with one of the two singular behaviours 
described in chapter H, i.e. those with a S-type singularity resp. with a pseudo- 
RN-type singularity; we denote them NAS resp. NAPRN solutions. We have 
determined several NAS families represented by the dashed-dotted curves of 
Fig. ^ The curve 1 staying to the right of the parabola coincides with the 



corresponding one found in p2[, whereas the others, staying essentially to the 



left of the parabola are new |2^. The basic NAS curve 1 intersects (once) 



only n = 1 (solid) curve for asymptotically fiat black holes. As will be ex- 
plained in chapter ||, the two NAS curves 6 and 7 accompanying the (dotted) 
NARN curve 5 are expected to merge with the NAS curve 3 close to = 0.9. 
Some of the NAS curves (e.g., 3 and 4) are expected to extend indefinitely 
to the right, but numerical difficulties (too violent "explosions") prevented 
us from continuing them further to larger values of Vh- They will intersect 
the (solid) curves for asymptotically fiat solutions with n = 2, 3, . . . zeros of 
W and therefore yield additional asymptotically fiat NAS black holes. 

Asymptotically (for big r/j) the basic NARN curve 2 approaches the basic 
NAS curve 1. Why this happens can be "understood" from Fig. |[ 

Finally there are the NAPRN solutions, which constitute a discrete set 



8 



according to the number of available free parameters at r = 0. We found 
several such solutions [23|. Few NAPRN solutions are shown in Fig. ^ and 



Fig. 1^. Only one of them has no maximum of r and was found in as 
well. Let us stress that although the NAPRN solutions do not correspond 
to asymptotically flat black holes they play essential role in the explanation 
on the cyclic behaviour of a generic solution in the EYM case (compare 
discussion in the next chapter |^) . 



5 Qualitative Discussion 



We shall now give a qualitative picture of the solutions and try to explain 
our numerical results. Let us briefly summarize what can be done (and what 
in fact has been done 



- get a qualitative understanding of the solutions 

- obtain a plateau - to - plateau formula in the EYM case relating quan- 
tities at one plateau (before "explosion" ) to those on the next plateau (after 
"explosion") pl 

- describe a simplified dynamical system, which reflects the main prop- 
erties of the generic solutions p2|, ^ 

Since the generic behaviour of the solutions is rather different in the cases 
with and without Higgs field, we shall treat the two cases separately. Let us 
first concentrate on the case without Higgs field. 

5.1 EYM theory 

For a "naive" understanding of the cyclic behaviour one can use a mechanical 
analogy. Introducing the "time" variable a = — ln(r) the EYM equations can 
be written in the form 

'"'^'"'--'^ -' "W, (13a) 





l + ^]B, (13b) 



where ' = d/dcr. The first equation Eq. ([T3|a)) resembles the motion of a 
ficticious particle in a potential with velocity dependent friction. 
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Note that the sign of the friction coefficient (term in square brackets in 
the Eq. (|l^a)) can be positive as well as negative. Close to the horizon 
(W^ — l)^/r^ — 1 < and friction coefficient is positive, corresponding to a 
deceleration of the "particle" . As the time a increases (r decreases) this term 
changes sign and the friction turns into anti-friction. The particle starts to 
accelerate quickly. This leads to a domination of the second term (kinetic 
energy) in Eq. (|T3|b), which in turn leads to a fast growth of the function B 
(respectively m). But growth of B stops the anti- friction in Eq. ([13|a) and 
the particle is again in the slow roll regime until the next "explosion" . 

For more detailed discussion of the generic behaviour we introduce the 
notation U = BW and B = rB and use again a = — ln(r) as a radial 
coordinate Note that B ^ —2m for small r. With these variables the 
field Eqs. (I) 



W = -r^^, (14a) 



B 

r B^ 



U = -W + 2r2^, (14b) 



B = r[^ ' ' -lj+2r^-g-. (14c) 

Close to the horizon the first term in the equation for B dominates (since 
f7 vanishes at r = r^) and thus B becomes negative. Provided W"^ does not 
tend to 1, this term will, however, change sign as r decreases and B will turn 
back to zero. Assuming further that IJ does not tend to zero simultaneously, 
the second term in the equation for IJ will grow very rapidly as B approaches 
zero, leading to a rapid increase of f7. This in turn induces a rapid growth of 
B (compare Fig. ^). Once the second terms in Eqs. (^^,c) dominate one gets 
{U / B) and thus IJ / B = W'/r tends to a constant c. As long as (rc)^ 
is sizable U and B increase exponentially, giving rise to the phenomenon 
of mass infiation. Eventually this growth comes to a stop when (rc)^ has 
become small enough. Then U and B stay constant until the first terms 
in Eqs. (p!^,c) become sizable again. As before B tends to zero inducing 
another "explosion" resp. cycle of mass infiation (compare Fig. H). 

In the discussion above we made two provisions - that stays away 
from 1 and that If does not tend zo zero simultaneously with B. If the 
first condition is violated, i.e. W"^ ^ 1 we get a NAS solution. If on the 
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other hand U and B develop a common zero we get a NARN solution, i.e. a 
solution with a second horizon. Both these phenomena can occur after any 
finite number of cycles, giving rise to several NAS resp. NARN curves as in 
Fig. 0. Generically W changes very little during an inflationary cycle, with 
the exception of solutions that come very close to a second horizon, i.e. close 
to a NARN solution. In this case W may change by any amount, depending 
on how small If becomes at the start of the explosion. By suitably flne-tuning 
the initial data at the horizon one can then obtain new NAS solutions with 
W ±1 or a new NARN solution. In this way each NARN solution is 
the 'parent' of two NAS and one NARN solution. This schematically shown 
on the Fig. ^. Fig. ^ shows two such generations: the NARN solutions 
labelled 2 have the NAS children 3 and 4 and the NARN child 5; the curves 
labelled 6 and 7 are the NAS children of 5. Whenever the value of W at the 
second horizon of a NARN solution approaches ±1 this NARN curve and its 
NAS children merge with the corresponding sibling NAS curve having one 
cycle less. This hierarchy of special solutions gives rise to a kind of chaotic 
structure in this region of "phase space" . 



Neglecting irrelevant terms one can integrate Eqs. and obtain a 



plateau - to - plateau relation |g3[, which connects the quantities Wo,Uo 
and Bq before an explosion with Wi, Ui, Bi after it 

[/i = t/oe(-o)\ 5i = ^e(-o)% W, = Wo-^rl, (15) 

with 

r„^JM^, c^^Mzll. (16) 



It is instructive to illustrate these relations on an example. We take the 
fundamental black hole solution with r/j = 1 and Wh = 0.6322 shown in 
Fig. 1 in P3|. For the flrst explosion one flnds the parameters ro ~ 2.7 ■ 10"'' 



and c ^ 1.1-10^ yielding ctq ^ 30 and thus Bi ~ e^°° and Wi — Wq^A- 10~^. 
The subsequent explosion will then take place at the fantastically small value 

Since the change of W in one inflationary cycle has an extra factor vq 
the function W stays practically constant. If we furthermore concentrate on 
cases, where the flrst term in Eq. (0b) can be neglected we may use the 

^Note that the similar relations are obtained recently in |H. 
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simplified system p2|, E3 



7-/3 . (I _ ]y2\2 tt2 

ij^ r B 

Introducing the variables x = rU / B = W and y = —(1 — l^^)^/ri? one 
obtains the autonomous system 

W = 0, x={y-l)x, y = y{y+l-2x^) . (18) 

Since the first of these equations may be ignored, we can concentrate on the 
x,y part. As usual for 2-dimensional dynamical systems the global behavior 
of the solutions can be analyzed determining its fixed points. Since the "large 
time" behavior a — oo corresponds to the limit r ^ these fixed points are 
related to the singular solutions at r = discussed in chapter ^ There are 
essentially three different fixed points. 

1. For y < there is the fixed point x = 0, y = —1 giving the RN 
type singularity. Its eigenvalues are —1 and —2, hence it acts as an 
attracting center for a —>■ oo. 

2. Then there is the point x = y = 0, a. saddle with eigenvalues ±1. 

3. In addition there are the points x = ±1, y = 1 with the eigenvalues 
1/2(1 ± iy/Tb), related to the pseudo-RN type singularity. This fixed 
point acts as a repulsive focal point, from which the trajectories spi- 
ral outwards. Since solutions of the approximate system given by the 
Eqs. (p!8D cannot cross the coordinate axes, solutions in the quadrants 
y > 0, X > resp. x < stay there performing larger and larger turns 
around the focal point coming closer and closer to the saddle point 



X = y = without ever meeting it. As observed in p2| this nicely ex- 
plains the cyclic inflationary behaviour of the solutions in the generic 
case. 

It is interesting to note that the similar results were obtained in the 
Abelian case p6| , p7[ in the homogeneous mass inflation model f\ 



^We are thankful to A.Ori for bringing the ref. to our attention and for communi- 
cating his unpubhshed results |M. 
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5.2 EYMH theory 

Finally let us discuss the black holes with Higgs field. Apart from the generic 
solutions there are the special ones approaching r = with a singular be- 
haviour described in chapter ^ On the other hand, the generic behaviour is 
much simpler than in the previously discussed situation without Higgs field. 
An easy way to understand this difference is to derive the analogue of the 
simplified system Eqs. ([18|) . Introducing the additional variable z = —H and 



ignoring again irrelevant terms one finds |23 



W = 0, H = -z, (19a) 
x = {y — l)x, z = yz (19b) 
y = y{y + l-2x'~z') . (19c) 

Leaving aside the decoupled equations for W and H one may study the fixed 
points of the {x, y, z) system. For z = one clearly finds the previous fixed 
points of the (x, y) system. However, for z 7^ the focal point disappears 
and the only fixed point for y > is x = y = z = zq with some constant 
zq. For < 1 this point is a saddle with one unstable mode, whereas for 
Zq > 1 it is a stable attractor. The latter describes the simple infiationary 
behavior described in chapter ^ and shown in the left part of Fig. ^. Solutions 
approaching a fixed point with Zq <1 eventually run away from it again and 
ultimately tend to one with > 1 as shown in the right part of Fig. ^. 
In the limit the equations can be trivially integrated with the result 



y = i/oe^^-^o^'^ = y^r^l~\ (20) 

Thus the mass function grows exponentially in "time" a or as a power in 
terms of r 

m = moe^o- = 1^ (21) 

in perfect agreement with our numerical results Fig. |^. Note that the fixed 
point Eq. (|20|) differs from the ones listed in the chapter ^ since the nature 
(strength) of the singularity depends on the solution itself. This fixed point 
was "rediscovered" in ||30[] . 

Note that the inclusion of a scalar field in the homogeneous mass infiation 
model 1^ has the same effect as the addition of the Higgs field to the EYM 



theory, namely that the mass infiation cycles disapear. 
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6 Concluding remarks 



To summarize the interior geometry of non-Abelian black holes exhibits a 
very interesting and complicated structure. Besides the generic solutions 
there are special NAS, NARN and NAPRN solutions, which can be obtained 
by a fine tuning initial data at the horizon. The main conclusion is that no 
inner (Cauchy) horizon is formed inside non-Abelian black holes in generic 
case, instead one obtains a kind of mass infiation. Without a Higgs field, i.e. 
for the EYM theory, this mass infiation repeats itself in cycles of ever more 
violent growth. 

A natural question to ask is what might be potential outcome of this 
investigation. A short answer would be 



illustration of singularity theorems |^ 



possible cosmological applications (see e.g. 



- interesting laboratory for non-perturbative study of (homogeneous) 
mass infiation phenomenon 

Naturally one should bare in mind limitations which come from quantum 
corrections and instability. Our consideration was purely classical, but when 
r tends to zero the curvature diverges and quantum corrections will become 
important. As it was discussed in the Introduction the EYM and most of 
the EYMH black holes are classically unstable. What will be a fate of a non- 
static perturbations inside non-Abelian black holes is an interesting open 
question. Definitely these subjects require further study. 
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Figure 1: Schematic view of the different special solutions in the EYM 
theory. 
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Figure 2: Initial data for special solutions. The solid curves represent 
asymptotically flat solutions with n zeros of W. The other curves represent 
various NARN and NAS families. 
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Figure 3: Parameters of the basic (lowest) NARN solutions. One clearly 
sees that with increasing rh the value of the second (inner) horizon rho — > 
and value of the gauge field there Who 1- This values correspond to NAS 
solutions. 
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Figure 4: First two cycles of the solution with r/i=0.97 and Wh = 0.2. For 
the second cycle a suitably stretched coordinate x is used. 
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Figure 5: Two characteristic types of inflationary solutions with Higgs fields. 
Note that in both cases asymptotically ln(m) is linear in ln(r). 
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Figure 6: Two different NAPRN solutions witli no zero of W between r — 
and r = r/j. 
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